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Abstract 

We revisit the gauge-covariant canonical formalism by separating explicitly physical and gauge degrees of freedom. 
We show in particular that the gauge-invariant linear and angular momentum operators proposed by Chen et al. can 
consistently be derived from the standard procedure based on the Noether's theorem. Finally, we demonstrate that this 
approach is essentially equivalent to the gauge-invariant canonical formalism based on the concept of Dirac variables. 
Because of many similarities with the background field method, the formalism developed here should also be relevant 
to general relativity and any metric theories. 



1. Introduction 

The canonical formalism allows one to derive the 
classical equations of motion and the conserved cur- 
rents associated with the global symmetries of a La- 
grangian. Unfortunately, this formalism runs into prob- 
lems in presence of a gauge symmetry. In particular, the 
canonical energy-momentum tensor obtained following 
the standard procedure turns out to be gauge dependent. 
Since physical observables are gauge invariant, there is 
a widespread belief that the canonical quantities are not 
really physical. Moreover, the gauge dependence of the 
canonical variables makes the quantization procedure 
particularly non-trivial. Similar problems arise in any 
metric theories like e.g. general relativity. 

Different strategies have been adopted to deal with 
these issues. The oldest one goes back to Dirac |l|] 
who proposed to reformulate QED in terms of gauge- 
invariant fields known as Dirac variables. These vari- 
ables are constructed by adjoining phase factors to 
the original fields, and have been rediscovered several 
times under different names and in different contexts 
ESI II ESS SHOP). We refer to Q for a 
review of the subject and to 1U4I1 for a comparison with 
the more standard Faddeev-Popov approach. 

Duan et al. EEI 171 11811 adopted a different strat- 
egy. They proposed to separate explicitly the physi- 
cal degrees of freedom from the unphysical ones in the 
gauge potential. The same idea reappeared in lfl9i Eoll 
and, more recently, in the context of the gauge-invariant 



decomposition of the proton spin 112 lL 1221 123L 1241 12511 . 
Interestingly, this approach shares many common fea- 
tures with the background field method introduced by 
DeWitt SHI SHI HI- The latter has been 
extensively used in gravity and supergravity 
361 37 , 3811. and in both conti nuum and lattice gauge 



theories SllSlinilllllSIIil], essentially in the path 
integral formulation. A nice introduction to the back 



ground field method can be found in 145 1. 

In this letter, we show how the explicit separation 
of physical and gauge degrees of freedom naturally 
leads to a covariant form of the Euler-Lagrange equa- 
tions and of the Noether's theorem, reconciling there- 
fore the canonical formalism with the gauge symme- 
try. In section [2] we briefly remind the textbook ap- 
proach to the canonical formalism. Then we present 
its covariant form in the presence of external or non- 
dynamical gauge fields, defining on the way covariant 
functional derivatives. In section [3] we develop the 
gauge-covariant canonical formalism based on the de- 
composition of the gauge field into physical and pure- 
gauge parts. Applying this new gauge-covariant canon- 
ical formalism to QCD, we recover the gauge-invariant 
decomposition of the linear and angular momentum op- 
erators constructed by Chen et al. 12 ll l22l 12411 . This 
naturally explains why their operators are the genera- 
tors of translations and Lorentz transformations for the 
physical fields. We also comment on the lack of unique- 
ness of this approach and its physical relevance. In sec- 
tion |U we discuss the Dirac's gauge-invariant canoni- 
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cal formalism, and demonstrate its equivalence with the 
gauge-covariant canonical formalism we propose. Fi- 
nally, we conclude this letter with section 

While we confine here our discussions to the gauge 
theories, the same approach can easily be adapted to 
general relativity and any metric theories. 

2. Canonical formalism 

We start with a short reminder of the standard deriva- 
tion of the Euler-Lagrange equations and the Noether's 
theorem. Then we explain how one reconciles the stan- 
dard approach with the gauge symmetry in presence of 
external gauge fields. 

2.1. Standard approach 

In standard textbooks like e.g. Hfl, the Lagrangian is 
usually thought of as a function of a generic set of fields 
4>{x) and their ordinary derivatives^] 



(1) 



Setting to zero the variation of the action 5S - 
J n d 4 x 6£(x) = under arbitrary (infinitesimal) varia- 
tion of the fields 



where £'(x') = f[<j>'{x'), d'^'ix')]. Using the Euler- 
Lagrange equations, one concludes that the following 
(infinitesimal) currents are conserved 



_d£_ 



A(f>- 



_d£_ 



8yd) - 5y£ 



5x v 



(6) 



When the Lagrangian is invariant under a gauge sym- 
metry, the Euler-Lagrange equations turn out to be 
gauge covariant. Note however that the terms gjr^j 

and |£ are not in general separately gauge covariant. 
More troublesome is the fact that the currents associated 
with the Poincare transformations are also not gauge in- 
variant. It is then often claimed that the canonical lin- 
ear and angular momentum operator densities have no 
physical significance. 



2.2. Gauge-covariant approach 

It is actually possible to reconcile the 



<p(x) M> (p'(x) - (p(x) + 6(f>(x) 



(2) 



auge symme- 
try with the Noether's theorem H l49l poi Pll pi P"3ll . 
The origin of the problem comes from the fact that 
the standard approach deals with quantities that are not 
gauge covariant, and therefore ill-defined from the geo- 
metrical point of view. 

To keep the presentation simple, consider that the 
fields (f>{x) transform as internal vectors under gauge 
transformations 



that satisfy the condition 6<p - on me space-time 
boundary dQ., one obtains the Euler-Lagrange equations 



d£ 
'3(cU 



d<p 



(3) 



In the Noether's theorem 04711 . one considers more gen- 
eral variations of the fields 



(x) ^ (p'(x') = (p(x) + A(j)(x) 



= <p(x) + 6<f>(x) + d^ix) dxf 1 , 



(4) 



where 6<p{x) represents again an intrinsic change in the 
functional form of the fields, d ll (p(x)6x fl represents a 
change coming from the fact that the fields are evalu- 
ated at a slightly displaced point x' — x + 6x, and Atp(x) 
is the total variation. By definition, continuous symme- 
tries leave the action invariant 



SS 



J 



d 4 x' £'(x) 



f d 4 x£(x) 
Jn 



0, 



(5) 



1 One can also add an explicit dependence on the space-time coor- 
dinates x, but this does not affect in a significant way our discussions. 

2 We adopt the passive point of view, so that the space-time volume 
is not affected. 



<p(x) H-» (f>(x) = U(x)(p(x). 



(7) 



All the following expressions can of course easily be 
adapted to any kinds of internal tensor transformation 
law. Since the original and transformed fields are eval- 
uated at the same point, it follows that the intrinsic vari- 
ation of the fields 6(p(x) transforms also as an internal 
vectoi[| under gauge transformations. The action being 
gauge invariant, one can deduce from 



6S 



-L 



d 4 x 



d£_ d£ 



(8) 



that the Euler-Lagrange equations are automatically 
gauge covariant. 

In a gauge theory, the ordinary partial derivatives are 
not the natural geometric objects but usually appear as 
part of covariant derivatives — - igA^. It is there- 
fore more natural to consider the Lagrangian as a func- 
tion of the fields and their covariant derivatives 



£{x) = f'[<j>{x),D^{x)l 



(9) 



3 Even when the transformation of the fields involves non-tensorial 
terms, the latter are cancelled in the intrinsic variation. 
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Simple algebra shows that |5£ 

df df df 



d<p d<p d(D^) 
df 8f 



did,,® d(D^y 



(10a) 
(10b) 



When the Lagrangian is expressed only in terms of 
gauge-covariant variables, the corresponding functional 
derivatives are automatically gauge covariant. It is 
therefore natural to define the covariant functional 
derivatives and conjugate fields as follows 



n',' = 



df 



df 



(lib) 



The Euler-Lagrange equations can then be rewritten as 



D u -— = 



8(D M <p) " dcf> 



(12) 



with the covariant derivative in the conjugate fundamen- 
tal representation given by D^-dfi +igA fl , or more com- 
pactly using the definitions (fTTT i 



(13) 



In this form, the individual terms of the Euler-Lagrange 
equations are now gauge covariant. 

The problem with the Noether's theorem in the stan- 
dard approach boils down to Eq. ©. When the fields 
transform e.g. as internal vectors under gauge transfor- 
mations, this expression does not make sense anymore 
from the geometrical point of view. Indeed, the fields 
</>'(x') and <p(x) live in different copies of the internal 
space since they are evaluated at different space-time 
points. It then follows that, contrary to the intrinsic vari- 
ation 6(f>(x), the total variation A<p(x) does not transform 
covariantly. The consistent expression is ||55, 561 



<p'(x') = W(x + 6x, x) [4>(x) + A c (f>(x)] , (14) 

where the infinitesimal Wilson line and the covariant to- 
tal variation are respectively defined as 

<W(x + 6x, x) = 1 + igA^x) 8xf, (15a) 
A c (f>(x) = 6(f>{x) + D^ix) fix?. (15b) 

From the gauge transformation of the A^(x) field 



A M (x) i— » A[i(x) = U(x) 



A^x) + -d ll 



U-\x), (16) 



it is easy to see that the infinitesimal Wilson line trans- 
forms in a simple way 

r W(x + Sx, x) i->'W / (x + 6x, x) 

= U(x + 5xYW{x + Sx, x)U~\x), (17) 

and allows one to parallel transport (p'(x') to <p'^{x) = 
<p(x) + A c <f>(x) that live in the same copy of the internal 
space as <p{x). The consistent expression for the con- 
served current is then 



_d£_ 

d(D4 



Ar <t> - 



d£ 



d(D^) 



D v d> - 



5x v . (18) 



This approach is fine as long as the gauge field is ex- 
ternal. But once A M (x) is treated as a dynamical field, 
the gauge-covariant formalism presented in this section 
cannot be applied anymore, simply because the gauge 
field does not transform as an internal tensor under 
gauge transformations. 

3. Gauge-covariant canonical formalism 

The gauge field plays essentially two roles. On the 
one hand, it is used to form a gauge-covariant derivative 
and to render the Lagrangian gauge invariant. On the 
other hand, it gives rise to a non-vanishing field strength 
= DfjAy - DyAft and provides the coupling with the 
source fields. The first aspect is somewhat unphysical 
as it concerns only the gauge symmetry which is not ob- 
servable. On the contrary, the second aspect is physical 
as the field strength (or curvature) affects the trajectories 
of particles, and is therefore observable. We adopt here 
the strategy of separating these two aspects explicitly. 

3.1. Decomposition of the gauge field 

Using the notation introduced by Chen et ai, we de- 
compose the gauge field as follows I2lll22ll24 l 



A„(x)=A*r(x)+A^Xx), 



phys/ 



(19) 



where A^ ure (x) and Ajf ys (x) contain only gauge and 
physical degrees of freedom, respectively. By defini- 
tion, the pure-gauge field is unphysical and therefore 
cannot contribute to the field strength 



ppuie 



d.Ar - d v A v ; K - ig [A£ ure , A?*] = 0. (20) 
It can then be written in the form 



A V ; a (x) = - f/pureM^t/pl^X (21) 
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where t/ pure (x) is some unitary matrix. From the gauge 
transformation law of this matrix 



(22) 



'purcW U puie (x) = U{x)U puK (x) 

and Eq. ( TTol l. it is easy to obtain the gauge transforma- 
tion laws of the pure-gauge and physical terms 



A P llre 



4 pure 



(x) 
U(x) 



U-\x), 



iphy; 



i phys/ 



(23) 
(24) 



In particular, note that A p (x) transforms as an internal 
tensor just like any other physical fields. It is therefore 
natural to treat the physical term as a dynamical field, 
i.e. as part the generic set of fields <p(x), and to treat the 
pure-gauge term as an external field. 

3.2. Gauge-covariant approach revisited 

Once more, consider for simplicity that the fields 
4>{x) transform as internal vectors under gauge transfor- 
mations. We propose to think of the Lagrangian as a 
function of the gauge-covariant physical fields and their 
pure-gauge covariant derivativesj D pme =8^- igA v ^ K 



£(x) = f"[<f>(x),DPr<Kx)l 



(25) 



Such a rewriting is always possible since a gauge- 
invariant Lagragian involves the field A /J (x) only 
through covariant derivatives and field strengths 

D^DVT-igAf*, (26) 
= Dl m ~ e A? hys - Dr e Af ys - ig [Af ys , A? hys ] , (27) 

where the pure-gauge covariant derivative in the adjoint 
representation is given by 2)JJ ure =8^- ig[A^ K , ]. 
Note that, owing to Eq. d20l i, the pure-gauge covariant 
derivatives commute. They are therefore the most nat- 
ural gauge-invariant extensions of the ordinary partial 
derivatives d^. 

In this approach the geometrically consistent expres- 
sion for the variation of the fields is obviously 

<P'{x') = <W VUK (x + Sx, x) [(f>(x) + A puK <p(x)] , (28) 

where the pure-gauge infinitesimal Wilson line and the 
pure-gauge covariant total variation are respectively de- 
fined as 



r W P ure(x + Sx, x) = 1 + igA p ™(x) 6x*, 



(29a) 



A VUK <p(x) = 6<p(x) + DjT>(*) &t". (29b) 



The Euler-Lagrange equations and conserved Noether 
currents then take the form 







d£ 



d£ 



■DL 



d£ 

d(f>' 



5x v 



(30) 



(31) 



d(D p ; ie 

Considering infinitesimal translations of the space 
time coordinates with the fields physically unchanged 



Sx? = ef, 



A pure 0W = 0, 



(32) 



we obtain from the corresponding (infinitesimal) 
Noether current = T^Sy the gauge-invariant canon- 
ical energy-momentum tensor 



jflV _ 



(33) 



Considering now infinitesimal Lorentz transformations 



with to 



ft)' 



field variations are given by 



A pure \ 



(x) 



the coordinate and 



j a pS a P<p(x), (34) 



where S a/3 is the appropriate (antisymmetric) spin ma- 
trix. From the corresponding (infinitesimal) Noether 
current J" = \ M^diyp, we obtain the gauge-invariant 
canonical generalized angular-momentum tensor 

MflVP = -'" ST^fc^ fl*) + (x v T» p - xPTn- (35) 

The crucial difference with the standard treatment is 
simply the use of the geometrically consistent total vari- 
ation A pure <P(x) instead of A(f>(x). 

There is therefore a simple rule of thumb for recon- 
ciling the standard canonical formalism with the gauge 
symmetry: it suffices to replace formally A^ by A phys 
and by the appropriate pure-gauge covariant deriva- 
tive in any gauge-dependent expression. 

3.3. Application to QCD 

Using the decomposition ( fT9l , the QCD Lagrangian 
can be thought of as made of three gauge-invariant terms 
-Cqcd = -Co + -Cym + -Cint, where the so-called Dirac, 
Yang-Mills and interaction terms are given by 



^pure(-r) = 0. 



X D = ^(iB pme - tri)tf/, 
X mt = g# phys <A- 



(36a) 
(36b) 
(36c) 
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We used for convenience the notation a = i (a - a). 

From the Euler-Lagrange equations (l30t , we recover 
the standard QCD equations of motion 



= 



tTpure 



o =d| 



pure 



34* 



= £>p 



dt 



d£ 



d(T% we Af ys ) dAf y " 



(37a) 



(37b) 



(37c) 



where a,b are internal-space indices. For the gauge- 
invariant canonical energy-momentum tensor (133K we 
obtain 

T HV = j^v + _ gfn ,^ (38) 

where the quark and gluon contributions are given by 



77 = itfri^, 

rf --2Tr[ J F^ ure AP hyS ]- 



(39a) 
(39b) 



Similarly, for the gauge-invariant canonical generalized 
angular-momentum tensor (l35l l. we obtain 



<y,spin <j,OAM 



+ m g,spin + iK, g,OAM 



where the spin and orbital angular momentum (OAM) 
contributions of quarks and gluons are given by 

KTpin = \ tWX^ = \ ^vysfc (4la) 



q,OAM 1 
g,spin L pnys J 

<OAM = - 2Tr [^ [V ^ureAP hyS ], 



(41b) 
(41c) 
(41d) 



with e ( )i23 = +1 and the notation a^b^ = a^V - a v lf. 

The expressions ( l39l and d4TT > coincide with the 
gauge-invariant canonical decompositions of the proton 
momentum and spin proposed originally by Chen et al. 
I 2ll 12211 . and put in a Lorentz-covariant form II24[ |25 1 . 
We therefore demonstrated here that the ad hoc expres- 
sions of Chen et al. can actually be derived from the 
canonical formalism and the Noether's theorem, once 
these are reconciled with the gauge symmetry. 



3.4. Stueckelberg symmetry 

By construction, the decomposition (TT9b is gauge in- 
variant, i.e. one has A^(x) = A^ ure (x) + AJ^ hys (jc). It is 
also consistent with the Lorentz symmetry, as discussed 
in detail in Ref. 12511 . However, it is not unique since 
we still have some freedom in defining exactly what 
we mean by 'pure-gauge' and 'physical'. The reason 
is that decomposing the gauge field into two parts au- 
tomatically introduces an additional local symmetry to 
the Lagrangian. The new symmetry has the same group 
structure as the gauge symmetry but acts only on the 
pure-gauge and physical parts of the gauge field 

AfTW + - U pwe (x)Uo\x) [d^Uoix)] U'Ux), 



(42a) 



Af ys ' s (x) 



Af y \x) - - U pm - e (x)U [ (x) [d M U (xj\ U p Ux), 

(42b) 

where Uq(x) is a gauge-invariant unitary matrix. At the 
level of the matrices i/ pure (x), this transformation reads 

f/ P ure(x) i-» U s pme (x) = U pme (x)U \x). (43) 

While the ordinary gauge transformation acts on the left 
of U puK (x) as in Eq. (l22l . this new transformation acts 
on the right. It is therefore important to distinguish 
them. Since the pure-gauge term AjJ ule (x) plays a role 
somewhat similar to the derivative of the Stueckelberg 
field, we refer to this transformation as the Stueckelberg 
(gauge) transformation 1E5I1 . Note also that the Stueck- 
elberg symmetry has no global counterpart, as one can 
see from Eq. (l42l . This means that the decomposition 
of the gauge field into pure-gauge and physical contri- 
butions does not introduce new conserved currents in 



the theory 157 1. 

The Stueckelberg symmetry is a bit problematic in 
the sense that one can write in principle (infinitely) 
many Lagrangians equivalent to the original one, just 
by changing the explicit expressions for A^ re (x) and 
A]J hys (x). To single out a particular Lagrangian in prac- 
tice, one can add a gauge-invariant term that breaks the 
Stueckelberg symmetry. This amounts to constraining 
further Ajf ys (x). One can use for example the light- 
front constraint A^ hys (x) = 0, the Coulomb constraint 
■gpuie . AW S ( X ) = 0, the Fock-Schwinger constraint 
x ■ A phys (x) = 0, or any other physical constraint that 
specifies the two physical degrees of freedom. It is im- 
portant to keep in mind that, despite appearances, this 
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procedure does not fix the gauge since A^ llle (x) also con- 
tributes to A M (x). 

Explicit realizations of the decomposition (fT9] l clar- 
ifies the physical meaning of the Stueckelberg symme- 
try. These realizations are usually non-local. The gauge 
invariance is then assured by the use of Wilson lines 
whose path dependence is at the origin of the Stueck- 
elberg symmetry J58|. In other words, breaking explic- 
itly the Stueckelberg symmetry amounts to determin- 
ing the path of the Wilson lines. This path dependence 
has a physical relevance as demonstrated e.g. by the 
Aharonov-Bohm effect lf59ll and the possibility to ac- 
cess the transverse canonical momentum of partons in 
the TMD factorization framework, see e.g. 16011 and ref- 
erences therein. 

It is the actual physical process that determines the 
shape of the Wilson lines, and in turn which constraint 
on A{f ys (x) to use. Wilson lines are necessary to pre- 
serve the gauge invariance, but their shape favors a 
particular gauge constraint, the one in which the Wil- 
son line factors reduce to the identity or, equivalently, 
the one in which A]J ure (x) vanishes. Any gauge will of 
course give the same numerical answer for the physical 
observable, but the physical interpretation of the latter 
will be the simplest in the gauge favored by the Wilson 
line. The archetypical example is deep-inelastic scatter- 
ing where the factorization theorem forces the Wilson 
lines entering the definition of the parton distribution 
functions to run along the light-front direction. At lead- 
ing twist, these parton distribution functions can be in- 
terpreted as linear combinations of parton probabilities 
in the light-front gauge A + (x) = 0. The decomposition 
( fT9b simply allows one to extend this interpretation to 
any gauge, provided that one defines the physical term 
by the constraint Ai (x) = 0. 



4. Gauge-invariant canonical formalism 

We present in this section an alternative to the gauge- 
covariant canonical formalism, which we refer to as 
the gauge-invariant canonical formalism. We show that 
these two formalisms are equivalent. 

4.1. Dirac variables 

Dirac soon realized that one of the main obstacles 
in the quantization of a gauge theory is the gauge de- 
pendence of the fields. He therefore built from the old 
gauge-variant fields new gauge-invariant fields that will 
play the role of the dynamical variables in the canonical 
formalism. In this spirit, we make use of the matrices 



U puK (x) to construct the field variables 
fax) = U^xMx). 



(44) 



For simplicity, we considered once more that the fields 
0(x) transform as internal vectors under gauge trans- 
formations, but this expression can easily be adapted 
to any sort of internal tensors. Similarly, for the fields 
like A jU (x) that do not transform as internal tensors under 
gauge transformations, we define 



AJx) = U'(x) 



A^x) + - d A 



(45) 



From the gauge transformation law (l22t of t/ pure (x), we 
see that the fields cf>{x) and A^(x) are by construction 
gauge invariant. 

We stress that, despite appearances, Eqs. (l44l 
and d45l > do not represent gauge transformations. In 
practice, the matrices I/ pure (x) can be expressed in terms 
of the gauge field A jU (x) 15811 . and can then be thought of 
as dressing fields. From a geometrical point of view, 
U vaK {x) simply determines a reference configuration in 
the internal space. The gauge-invariant fields 0(x) then 
represent "physical" deviations from this reference con- 
figuration. 

4.2. Gauge-invariant approach 

In the gauge-invariant canonical formalism, the La- 
grangian is thought of as a function of the Dirac vari- 
ables 4>{x) (including from now on A jU (x) as well) and 
their ordinary derivatives 



£(x)=/"'[0(x),<9^(x)]. 



(46) 



Since the Dirac variables are gauge invariant, we can 
simply apply the standard approach of Section |2TI 

The rule of thumb is particularly simple: it suf- 
fices to add a hat to the fields <p(x) wherever they ap- 
pear. The obtained Euler-Lagrange equations and con- 
served Noether currents are then automatically gauge- 
invariant. 

4.3. Equivalence with the gauge-covariant approach 

The gauge-covariant and invariant canonical for- 
malisms are essentially equivalent. Indeed, noting thajf] 



8^(x) = f/ pu 1 re (x)DP ur V(x), 



(47) 



For the gauge field, we have A^x) = U p J„,(x)Af, ys '(x)U pm .(x) 
-1 

pure 



anddwV*) = u£ a (x)[!l% m Af y, (xj]Up a *(x). 
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the Lagrangian can be rewritten in the following form 
£(x) = /"' [f/ pu 1 re (x)0(x), U-Ux)D*r<f>(x)~] ■ (48) 

Then, thanks to the gauge symmetry of the Lagrangian, 
we are assured that all the matrices U pme (x) disappear in 
the final expression so that we can write the Lagrangian 
in the form of Eq. d25l l. 

In other words, one can switch between gauge- 
covariant and invariant canonical formalisms by a mere 
change of variables. This also means that the issue of 
uniqueness raised by the Stueckelberg symmetry affects 
the gauge-invariant canonical formalism as well. In- 
deed, we see from Eq. d43l l that, even if the Dirac vari- 
ables are gauge invariant, they are not Stueckelberg in- 
variant 

(fix) i-» = l/otote (49) 

There is the same freedom in defining precisely 4>{x) as 
in defining A^, hys (x). The existence of entire class of 
composite fields was already pointed out by Dirac and 
Steinmann EES]. 

In particular, we note that the Dirac's gauge-invariant 
formulation of QED is equivalent to the Chen et al.'s 
gauge-covariant formulation, since both make use of the 

explicit dressing matrices U pu[e (x) — e ^- leading to 
the Stueckelberg-fixing constraint V ■ A phys (x) = 0. This 
particular choice makes the Coulomb gauge V -A(x) = 
special. In that gauge, the gauge-fixed fields coincide 
with the gauge-invariant ones 0(x)|^ a, _ = (p(x). For 
this reason, the gauge-covariant and invariant canonical 
formalisms can be interpreted as gauge-invariant exten- 
sions of the standard canonical formalism. 

5. Conclusion 

In this letter, we showed that separating explicitly 
the physical and unphysical degrees of freedom in the 
gauge potential allows one to reconcile in a natural 
way the Euler-Lagrange equations and the Noether's 
theorem of the standard canonical formalism with the 
gauge symmetry. Applying this formalism to QCD, we 
derived canonically the gauge-invariant operators pro- 
posed earlier by Chen et al. in the context of the pro- 
ton spin decomposition. Finally, we demonstrated the 
equivalence between the our formalism and the Dirac's 
gauge-invariant canonical formalism. 

Because of the similarity between the approach 
adopted here and the background field method, we be- 
lieve that the formalism developed in this letter should 
also be relevant to general relativity and, more gener- 
ally, to any metric theories. 
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